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High auditory sensitivity, sharp frequency selectivity, and otoacoustic emissions are signatures 
of active amplification of the cochlea. The human ear can also detect very large amplitude sound 
without being damaged as long as the exposed time is not too long. The outer hair cells are 
believed as the best candidate for the active force generator of the mammalian cochlea. In this 
paper, we propose a new model for the basilar membrane oscillation which successfully describes 
both the active and the protective mechanisms by employing an energy depot concept and a critical 
velocity of the basilar membrane. One of the main results is that thermal noise in the absence 
of external stimulation can be amplified leading to the spontaneous basilar membrane oscillation. 
The compressive response of the basilar membrane at the characteristic frequency and the dynamic 
response to the stimulation are consistent with the experimental results as expected. Our model 
also shows the nonlinear distortion of the response of the basilar membrane. 

PACS numbers: 05.45.-a, 87. 19. In, 87.16.Uv 



I. INTRODUCTION 



Remarkable detection capabilities appearing in animal hearing are essentially governed by both a passive mechanical 
and an active biophysical procedures in the cochlea [l|, 0; S 13 ■ The active amplification of the living cochlea was 
conjectured by Gold in 1948 Q and is now qualitatively studied widely through the Hopf equation in mathematical 
models 0, 0, M^ B [Tot - In order to give a better understanding to the experimental observations, a more accurate 
model is required. Furthermore the physical origin of the active oscillation remains still unclear. In mammals, the 
outer hair cells (OHCs) are known to be force generators for auditory sensitivit y a nd frequency selectivity. Proposed 
mechanisms of the force generation are contraction of the OHC itself [ll|, [l^, Tisj and an active motion of the hair 
bundle ^M, 15]. 

Energy depot model can describe the active phenomena [H, [l^, because the energy supplied by the depot can 
induce a neg ative dissipation. A similar mechanism for a negative stiffness in the bullfrog's hair bundle has been 
reported [131 ■ It is reasonable to consider the OHC as an energy depot since the OHCs play an important role in 
the active amplification of the cochlea even though we do not know the exact mechanism of the energy supply. In 
nature all living things have some mechanisms of protection for survival. We expect that there might be a similar 
mechanism to protect cochlea from damage by any external factors such as loud stimulation. Hence, we propose a 
braking mechanism in the cochlea to prevent the basilar membrane (BM) from an excessive oscillation which could 
damage the cochlea. For this purpose, a critical velocity of the OHC oscillation will be introduced in our model. 

In the following sections, we describe our model and reproduce the essential phenomena of the active cochlea 
observed in recent experiments such as the compressive nonlinearity, the dynamics of the BM response, nonlinear 
distortion and the spontaneous BM oscillation (SBMO). It will also be pointed ut that our model naturally includes 
the Hopf bifurcation model in the case of weak stimulation. However, in the regime of strong stimulation, the 
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present model is essentially different from the Hopf bifurcation model. 
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FIG. 1: The compressive response of the BM. The CF is lOkHz. It shows that the response of the BM is non-linear and 
compressive at the CF. (Note that the slope of the solid line is less than one.) Otherwise, the response of the BM is a linear 
function of the stimulation yielding that the slopes are 1. These results agree with the experimental observations [24| . 



II. ENERGY DEPOT MODEL WITH BRAKING MECHANISM 



It is well known that a change of transmembrane potential in the OHC induces a somatic vibration. This implies 
that a part of the electric energy is converted into the mechanical energy. In this sense, we believe that the OHC 
functions as an energy depot for the BM. The OHC as the energy depot stores energy supplied and converts it into 
kinetic energy overcoming dissipation. This energetics can be described by the following energy balance equation . 

'^=q-cE{t)-d{v)E{t), (1) 



where E{t) is the energy density of the energy depot, q the rate of the energy pumping into the energy depot per 
unit area, c the rate of the internal energy dissipation, d(v) the rate of the energy converting into the kinetic energy 
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of the BM, and v the velocity of the BM. Considering that constant contribution of d{y) can be incooperated into 
c and that conversion rate of the energy into the kinetic energy of the BM wiU not be sensitive to the sign of w, we 
assume that d{v) is an even function of v without a constant term. The energy depot model in which d{v) = d2v'^ 
has been extensively discussed fl6j . In this paper, we include the next higher contribution to d{v) and introduce a 
critical velocity of the BM to describe a braking mechanism [17.]. Hence the conversion rate d{v) can be written as 

d{v) = d2V^ (l - 5) ' (2) 

where Vc is the critical velocity. For a finite Vc the oscillation of the BM can be braked by the negative conversion 
rate into the kinetic motion of the BM. Note that Vc ^ oo corresponds to a system without the braking mechanism. 
The critical velocity will be roughly estimated from the contraction of the OHC later for our numerical calculation. 
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FIG. 2: The dynamic response of the BM when the CF is lOkHz and the stimulation is lOOdB. The onset time is much 
shorter than the offset time which is not shown in full scale. This result is qualitatively consistent with the experimental 
observations /i — 600kgm~^s~^, d2/c = 6 x 10^m~^s^, and qd2/c — /i. The stimulation is unloaded at t = 9ms. The 

lower panel shows the response represents the case without the braking mechanism, showing a continuous increase as long as 
the stimulation is on. 



Now it is worthwhile to compare our newly proposed model with the one previously introduced in Ref. [lO|. 
Whereas we treat the energy density {E(t)) of the energy depot(OHC), in Ref. 10] the energy density (e(t)) of the 
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FIG. 3: (color online) The dependence of the response of the BM on four different sets of parameters, /i = GOOkgm ^, 
d2/c = 6 X lOV'^s^, Vc = O.OOlm/s (solid thin line); ^ = 6, 000kgni"^s"\ d2/c = 6, 000m"^s^ Vc = O.Olm/s (solid thick line); 
/i — 600kgm^^s~^, d2/c = 6 x lO^m^^s^, Vc — oo (dashed thin line), and fi = 6, OOOkgm^^s"^, d2/c = 6,000m~^s^, Vc — ao 
(dashed thick line). The velocity of the BM is (not) saturated when the stimulation is sufficiently large with (without) the 
braking mechanism. 



propagating hydrodynamic wave is considered. The model in Ref. [T3] describes the coupling of the active elements 
to the propagating wave, while the present model describes the coupling between the energy depot (OHC) and the 
BM oscillation. The emphasis in this paper is on showing that the active and compressive response arises naturally 
from the generalized energy depot concept. 

III. RESULTS AND DISCUSSION 

An adiabatic approximation in which the adaptation of the energy depot is very fast [l6j yields that 

= n 2/ 2^ 2 - (3) 

C + 0,2 (1 — /v^) v' 
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It should bear in mind that E{t) is not constant. The osciUation of the BM is then governed by 

pv + p.V + kR ^ Fpassive + FacUve, (4) 

where p is the density, p the physical damping coefficient, k the stiffness, and R the displacement of the BM. Fpassive 
includes the contributions of a sound wave (external stimulation) and a noise. The noise is assumed to be a Gaussian 
white noise. FacUve is the active force density acting on the BM by the OHC, 

_ d{v)E{t) 'c^v 

^active — , / 9\ ' \^ ) 

V do I , V 



i>2 



Note that the Langevin equation, Eq. describes the passive and the active response in a unified and natural way. 
It should be noted that the active force, Eq. ([5]), can be either positive or negative depending on the value of v. 
When is smaller than (weak stimulation) , a positive active force is provided on the BM (active mechanism) . On 
the other hand, when is larger than (strong stimulation), a negative active force, hence an extra drag force, is 
exerted on the BM (braking mechanism). 

The mass density of the BM is known to be about 0.77kg/m^ [T9| and the reported damping coefficient, p, of the 
BM has a large variation from 600 to 63, 000kg m^^s^^ [l^, llo] depending on the models. For a numerical calculation, 
we set the damping in the range of 600 6,000kg m~^s~^. It has been reported that the Guinea pig OHCs are able 
to change its length up to 5% when the transmembrane potential is varied. This corresponds to a displacement of 
1/im for the OHCs of 20/im at the high-frequency end [2l!]. Hence, in the regime of high frequency of lOkHz, the 
adjustable parameter Vc is roughly estimated about O.Olm/s. By considering the OHC as a capacitor and using the 
observed membrane potential, membrane capacitance [l^, and the scale of the OHC [l^, we estimate the value of 
q to be about O.Olkg s~^. The response of the BM to the sound wave can be then obtained straightforwardly using 
Eqs. ([U and ([5]). In the following calculations, we set qd^jc = p and the characteristic frequency (CF), — \f^Jp- It 
will be shown later that these relations arise from the conditions of the Hopf bifurcation. In the present calculations, 
the noise width of the distribution is chosen to be zero dB. 

The compressive and dynamical responses of the BM are obtained as shown in Fig. [T] and Fig. [21 respectively. Fig.[T] 
demonstrates that the response of the BM is non-linear and compressive at the CF, while the response is linear at other 
frequencies. This compressive response of the BM agrees well with the observations [1, as expected. Fig. [5] shows 
the time course of the BM response, reproducing the experimental observation that the onset time of the response 
is shorter than the offset time, although the onset time is somewhat larger than the experimental values [23|. It is 
important to notice that the velocity of the BM increases continuously as long as the stimulation is being turned on if 
the braking mechanism is not included (vc — oo). Hence, the braking mechanism introduced in our model is essential 
to explain the experimental observation at a sufficiently strong stimulation. A dependence of the BM velocity on 
Vc is plotted in Fig. [31 The results clearly show that lower critical velocity produces smaller response at the same 
stimulation, thus producing a supercompressive behavior as observed in experiment [l^]. This is the effect of the 
braking mechanism, which is significantly different from the Hopf bifurcation model. Note that the saturation of the 
BM response disappears if the braking mechanism is not considered. 

The gain of the oscillator is defined as the ratio of the displacement of the BM to the stimulation, 

^passive 

Using Eq. ((3|), we obtain the gain, which is a scaling function of the stimulation when the stimulation is not too 
large, G ~ P^ast^ie- This result agrees well with the experiment and the result of the Hopf bifurcation model 0], 

G ^ Fpassvia- Indeed our model includes the Hopf model naturally. To show this, we expand the active force, Eq. (O, 
up to the lowest nonlinear term when ^2(1 — /vDv^ /c is small. Eq. ^ is then approximately written as 

Fpassive = PV + (^P - q^^ V + kR 

d2 f I d2\ 2, /r,^ 
+ 9— - + - P ■ 7 
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FIG. 4: The amplification of tlie tfiermal noise in the absence of external stimulation. When qd2 /c < ^, the Brownian fluctuation 
is suppressed, while in the region of qdijc > /i, the Brownian fluctuation is amplified. For the plots, the intensity of the 
thermal noise is OPa, and the strength of the fluctuation is OdB. (a) ^ = 6, 000kgm~^s~^, Vc = O.Olm/s. (b) [i, — 600kgm~'^s~^, 
= O.Olm/s. 

In the Fourier spectra space, this directly corresponds to the Hopf equation [6*] with the following bifurcation conditions, 

1^- (9) 

c 

Hence our energy depot model with a braking mechanism reduces to the Hopf bifurcation model in the limit of weak 
stimulation. 

One of the most important results of our model is an amplification of the thermal noise itself in the absence of an 
external stimulation. At the bifurcation point or when qd^jc < the thermal noise is well suppressed. However, 
when qd2/c > /i, so that too much of the internal energy is converted into the kinetic energy, the thermal noise can 
be largely amplified as shown in Fig. [4l Here, we set the mean average of the noise at OPa and the width at OdB. 
Since it can be generally assumed that the oscillators are distributed in the vicinity of the bifurcation point, the 
noise may be amplified incoherently in the absence of the input signal. However, when a weak input signal with the 
same CF is introduced, it may phase-lock to the already existing amplified noise, thus enhancing the selectivity as 
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FIG. 5: Nonlinear distortion of the response of the BM. /i — 7.89kHz, /2 — 8.28kHz, both strengths of the stimulation are 
lOOdB, the response is measured at the place where the CF = 7.89kHz, — 10~^m^/s'^ and fi — 600kgm~'^s~^. 



discussed in Ref. [25| . Because it is believed that the SBMO is crucial for understanding the spontaneous otoacoustic 
emission [2^, our model may provide a clue to explore this interesting and important phenomenon. Fig. O shows the 
nonlinear distortion of the response of the BM at the position where the CF is 7.89kHz. The stimulation contains 
two tunes, /i — 7.89kHz and /2 — 8.28kHz. The strengths both are lOOdB. Our model shows clearly the distortions 
at 2/i — /2, 2/2 — /i in the response of the BM. 



IV. SUMMARY 



In this paper, we propose a theoretical scheme for the active and passive response of the BM based on a concept of 
the energy depot model with a braking mechanism. According to the experimental observations, the OHC is assumed 
to play a role of the energy depot by supplying an active force to the BM and also by reabsorbing an excessive kinetic 
energy like the regenerative-brake system in automobiles and electric vehicles [l^l • It is shown that a single equation 
of motion, Eq ([5]), can produce all the essential passive and active behaviors. Moreover, our model shows that the 
Brownian noise can be sufficiently amplified under a certain condition, thus leading to the SBMO. Although we applied 
our model only to the mammalian cochlea, it is equally possible to be applied to the nonmammalian vertebrates [2^ 
because the only required physical characteristics are a supply of the active force and a self-adaptation. In this sense, 
we believe that the concept of the energy depot model with a braking mechanism is quite universal and can be applied 
to any underdamped biological systems. 
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